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(STAR Collaboration) 

We present high precision measurements of elliptic flow near midrapidity (|i/| < 1.0) for multi- 
strange hadrons and f) meson as a function of centrality and transverse momentum in Au-|-Au 
collisions at center of mass energy ^sjvjv = 200 GeV. We observe that the transverse momentum 
dependence of 0 and Q. V 2 is similar to that of vr and p, respectively, which may indicate that 
the heavier strange quark flows as strongly as the lighter up and down quarks. This observation 
constitutes a clear piece of evidence for the development of partonic collectivity in heavy-ion collisions 
at the top RHIC energy. Number of constituent quark scaling is found to hold within statistical 
uncertainty for both 0-30% and 30-80% collision centrality. There is an indication of the breakdown 
of previously observed mass ordering between 0 and proton V 2 at low transverse momentum in the 
0-30% centrality range, possibly indicating late hadronic interactions affecting the proton V 2 - 


At sufficiently high temperature and/or high density 
Quantum Chromodynamics (QCD) predicts a transition 
form hadronic matter to de-confined partonic matter [l| . 
The main goal of the STAR (Solenoid Tracker at RHIC) 
experiment at the Relativistic Heavy Ion Collider 
(RHIC) is to study the properties of QCD matter at 


extremely high energy and parton densities, created 
in the heavy-ion collision. In high energy heavy-ion 
collisions, particles are produced with an azimuthally 
anisotropic momentum distribution, which is a result of 
hydrodynamical flow of the Quark-Cluon-Plasma (in the 
soft regime). One way to examine this anisotropy is to 
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measure elliptic flow (U 2 ), which plays a crucial role in 
the study of the QCD matter formed during the collision. 
The elliptic flow, defined as V 2 = (cos2((/? — 'I')), is the 
second Fourier coefficient of the azimuthal distribution 
of the emitted particle with respect to the reaction plane 
(defined by the beam axis and a vector between the 
centers of the colliding ions). Here (p is the azimuthal 
angle of emitted particle and 'I' is the azimuthal angle of 
the reaction plane. Over the past decade, experimental 
measurements have shown elliptic flow to be especially 
sensitive to the initial phase and equation of state 
of the system formed in heavy-ion collisions 
However, information about the early dynamics of the 
system may be modified by hadronic re-scattering in 
the later stage of the collision [3, |^. The hadronic 
interaction cross-sections of <(), S and H are expected to 
be small Q and their freeze-out temperatures are close 
to the quark-hadron transition temperature predicted by 
lattice QCD ©[HI- Hence, these hadrons are expected 
to provide information primarily from the partonic stage 
of the collision [l2Ulfl . Previous measurements of 4> and 
H V 2 from STAR [13] were statistically limited and little 
is known about the centrality dependence of H V 2 - The 
measurements of (j) and D V 2 presented here as a function 
of both transverse momentum (pr) and centrality help 
to alleviate these limitations. Moreover, high precision 
measurements of ())-meson V 2 relative to proton V 2 at low 
Pt may provide information on the effect of hadronic 
re-scattering in the late stages of the collision. 

We present the collision centrality and pr dependence 
of the elliptic flow of 7 r+ -|- 7 r“, K'^ + K~, Kg, p + p, 

(j), A -I- A, -I- s'*" and -|- D"*". For this study we 
used 730 million of Au-|-Au events at ^snn = 200 
GeV recorded by STAR in 2010 and 2011 with a 
minimum-bias trigger [l8l| . The collision centrality is 
determined by comparing the measured raw charged 
hadron multiplicity from the Time Projection Chamber 
(TPC) within a pseudorapidity window I 77 I < 0.5 with 
Glauber Monte-Carlo simulations [l^, The TPC 

and Time of Flight (TOF) detectors with full azimuthal 
coverage are used for particle identification in the central 
rapidity region (jr;] <1.0 for TPC and jTyj < 0.9 for 
TOF). Charged particles are identified using specific 
ionization energy loss as a function of momentum (in 
the TPC) and square of the particle mass as a function 
of momentum (for the TOF). We reconstruct short-lived 
iCg, A, S, H and (f) through the following decay channels 
: iCg —>'7r+-|-7r“,A—j-p-l-Tr, S—^-A -I-tt, Q ^ A + K 
and (j) —>■ A+ + K~. Topological and kinematic cuts are 
applied to reduce the combinatorial background for K'g, 
A, S and H. The detailed description of the procedures 
can be found in Refs. [illlia 
The rj sub-event plane method 13 is used for the elliptic 
flow analysis. An 77 gap of jAr;] > 0.1 between positive 
and negative pseudorapidity sub-events is introduced to 
suppress non-flow effects. The V 2 for short-lived hadrons 
(ATg, (j), A, S and H) is calculated as a function of 
invariant mass for each px and centrality bin in order 


to take into account the invariant mass dependence of 
the signal to background ratio. Details of this method 
can be found in Ref. [l^l- The observed V 2 values 
are corrected for finite event plane resolution which is 
determined by comparing the two 77 -sub event plane 
angles. A resolution correction is done by dividing the 
term cos 2 (<77 — '£' 2 ) by the event plane resolution for 
the corresponding centrality for each event following the 
method described in Refs. [13, [13 ■ The change in V 2 
between present method of resolution correction and the 
previous method used in earlier STAR publication [iSj 
is < 5% at y/sNN = 200 GeV. Here 'I '2 is the 2^^^^ order 
event plane which is used for V 2 measurements. 

For all particle species, the cuts used for particle 
identification (PID) and background subtraction are 
varied to estimate the systematic errors. Furthermore, 
different techniques (e.g. by counting entries in each bin 
of the invariant mass histogram or by fitting the shape 
of the invariant mass distribution using a function) 
for yield extraction are used. For tt^, and p{p), 
6 different combinations of track cuts and 3 different 
sets of PID cuts which finally yield 18 combinations 
have been used. For other strange hadrons {Kg, (j). A, 
S and H) 20 different cut combination are used. The 
root-mean-square value of point-by-point difference from 
the default value {v 2 from default set of cuts) is used 
as the systematic error on each data point. The total 
systematic error depends on px, centrality and particle 
species. We observed 3-5% systematic error for px < 1.5 
GeV/c and 0-30% centrality for (j), Kg, A whereas for 
S and D the systematic error varies from 8 % to 14%. 
Total systematic errors are less than 1% for tt^, K^ and 
p{p) for all Px and centralities. 

We investigated the effect of track reconstruction 
efficiency on the measured V 2 of identified hadrons 
for wide centrality bins, such as a 0-80% centrality 
bin. The centrality dependence of track reconstruction 
efficiency biases the measured V 2 toward events with 
higher reconstruction efficiency, an effect we will refer 
to as an “efficiency bias”. Due to the efficiency bias, 
the V 2 of S and D, each having three daughters, changes 
by no more than 5% in 0-80% centrality. For the other 
measured particles, the effect is less than 3% for 0-80% 
centrality. The V 2 {px) of all particles presented here 
have been corrected for the efficiency bias by using the 
inverse of efficiency as a weight for the V 2 as a function 
of Px and centrality. 

An additional correction is needed for (j), S and D 
V 2 - An event bias is naturally introduced when one 
measures V 2 in wide centrality bins, especially for 
the rare particles. As the measured V 2 is an average 
over all events weighted by particle yield, the average 
event shape depends on the particle type. A Glauber 
model [13 study of the average initial participant ec¬ 
centricity indicates the multi-strange hadron V 2 is more 
biased toward central events than that of the light and 
strange hadrons. Specifically, the average eccentricity 
for multi-strange hadrons in wide centrality bins is 
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FIG. 1. (Color online) The V 2 as a function of pt near midra- 
pidity (|i/| < 1.0) for (a) E~ + H"*" (b) Q~ + and (c) (j) 
from An + An collisions at y/SNN = 200 GeV for 0-30% and 
30-80% centrality. The systematic uncertainties are shown 
by shaded boxes and the statistical uncertainties by vertical 
lines. 


smaller than the eccentricity determined by the particle 
yield of all charged hadrons. One should take this effect 
into consideration if any conclusion on the number of 
constituent quark scaling is drawn. This bias can be 
corrected by normalizing the measured V 2 by the ratio of 
eccentricity to that weighted by the yield of the particle 
of interest. We find the event bias correction factors for 
0-30%, 30-80%, and 0-80% centralities are 1.002, 1.053 
and 1.028 for cj); 1.019, 1.054 and 1.091 for S; 1.068, 
1.067 and 1.177 for ft. The event bias correction for 
light and strange hadrons is small (< 0.03), perhaps 
be due to their copious production. Therefore, in the 
later discussion of number of constituent quark (NCQ) 
scaling, the event bias correction is applied only to the 
V 2 of multi-strange hadrons and (j) meson. The above 
correction factors remain almost unchanged if we use 
Color-Glass Condensate (CGC) based model to 
calculate eccentricity. 

In Figure [T] we present the elliptic flow parameter 
V 2 {pt) at midrapidity {\y\ < 1.0) for (a) + S'*", (b) 

+ O and (c) (j) in An -|- An collisions at ^/sf^ = 
200 GeV for 0-30% and 30-80% centrality. Event bias 
correction factors have been applied to the results shown 
in Fig [T] A clear centrality dependence of V 2 {pt) is 
observed for (p, S and fl, similar to that of identified light 
and strange hadrons previously measured by the STAR 
experiment [1^. The values of V 2 are found to be larger 
in peripheral collisions (30-80% centrality) compared 
to those in central collisions (0-30% centrality). This 
observation is consistent with an interpretation in which 
the hnal momentum anisotropy is driven by the initial 
spatial anisotropy. 

The NCQ scaling in V 2 for different identified hadrons is 
considered to be a good probe for studying the strongly 
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FIG. 2. (Color online) The V 2 as function of pr for tt, p (panel 
a) and cj), Q (panel b) from minimum bias Au-|-Au collisions 
at ^/sNN = 200 GeV for 0-80% centrality. The systematic 
uncertainties are shown by the shaded boxes while vertical 
lines represent the statistical uncertainties. 


interacting partonic matter. The observed NCQ scaling 
of identified hadrons in experimental data indicates 
the importance of parton recombination in forming 
hadrons in the intermediate pr range ( 2.0 GeV/c < 
Pt < 4.0 GeV/c) [ 30 l - [^ . Such scaling may indicate 
that collective elliptic flow develops during the partonic 
phase. Previous measurements have found that V 2 of tt, 
K, p, Kg, A, S and (p follow NCQ sca ling well at top 
the RHIC energy {■,,/snn = 200 GeV) [ 2 ^. The large 
statistics data sets collected by STAR detectors allow 
us to measure elliptic flow of multi-strange hadrons, 
specifically that of the baryon which is made of pure 
strange (s) or anti-strange (s) constituent quarks and of 
the p meson, consisting of one s and one s constituent 
quark. 

Figure [ 2 ] shows the V 2 as a function of pr for tt, p, 
p and n for 0 - 80 % centrality in An - I - Au collisions at 
y/sNN = 200 GeV. Here p and V 2 are corrected for 
the event bias mentioned earlier. Panel (a) of Figure [2] 
shows a comparison between V 2 of tt and p, consisting 
of up (u) and down (d) light quarks, and panel (b) 
shows a comparison of V 2 of p and H containing heavier 
s quarks. The V 2 of p and H are mass ordered at 
low Pt and a baryon-meson separation is observed at 
intermediate pt- It is clear from Figure [2] that the 
V2{pt) of hadrons consisting only of strange quarks {p 
and n) is similar to that of tt and p. However, unlike tt 
and p, the p and H do not participate strongly in the 
hadronic interactions, which suggests that the major 
part of collectivity is developed during the partonic 
phase in Au + Au collisions at y/s^N = 200 GeV. 

We now compare our results for NGQ scaling for differ¬ 
ent collision centrality classes to see how the partonic 
collectivity changes with system size. Figure [3] shows 
the V 2 scaled by number of constituent quarks (riq) as 
a function of pT/ng and {rriT — mo)/nq for identified 
hadrons from Au + Au collisions at ^/snn = 200 
GeV for 0 - 30 % and 30 - 80 % centrality, where rriT and 
niQ are the transverse mass and rest mass of hadron, 
respectively. Here, p, S and H V 2 are corrected for the 
event bias mentioned above. To quantify the deviation 
from NGQ scaling, we fit the Kg V 2 with a third-order 
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FIG. 3. (Color online) The 1)2 scaled by number of constituent quarks {uq) as a function of and {mT — mo)/nq for identified 

hadrons from An + Au collisions at ^satjv = 200 GeV. Ratios with respect to a polynomial fit to the Kg V 2 are shown in the 
corresponding lower panels. Vertical lines are statistical uncertainties and shaded boxes are systematic uncertainties. 


polynomial function. We then take the ratio of V 2 for 
the other measured hadrons to the Kg fit. The ratios 
are shown in the lower panels of Figure [3] Table |T] shows 
the deviations of cj), A, S and 52 V 2 from the Kg ht line 
in the range (mr — mo)lnq > 0.8 GeV/c^. 

For both 0-30% and 30-80% centralities, the scaling 


TABLE I. Deviation from the Kg fit line in the range {rriT — 
mo)fnq > 0.8 GeY/c? for 0-30% and 30-80% centrality. 


Particle 

Deviation 

0-30% centrality 

30-80% centrality 

b 

2.7±2.6(stat.)±1.8(sys.)% 

1.2±1.3(stat.)±0.6(sys.)% 

A 

4.3±0.8(stat.)±0.2(sys.)% 

1.5±0.7(stat.)±0.2(sys.)% 

H 

11.3±2.3(stat.)±1.4(sys.)% 

8.5±2.0(stat.)±0.5(sys.)% 


10.1±8.4(stat.)±5.3(sys.)% 

7.0±6.0(stat.)±1.5(sys.)% 


holds approximately within 10%, excluding pions. The 
deviation of pions could be due the effect of resonance 
decay and non-flow correlations [s^. We have seen 
similar order (^10%) of deviation if we use pr/nq 
scaling as a reference. The maximum deviation from 
NCQ scaling is ^20% at /y m = 2.76 TeV as observed 

by ALICE experiment [^. Therefore, at top RHIC 
energy, NCQ scaling holds better than LHC energy. The 
observed difference between the charged kaon and Kg V 2 
at low pt is due to differences in the pile-up protection 
conditions used in collecting the different data sets. The 
difference is taken to be an additional contribution to 
the systematic error on Kg V 2 - 
Hydrodynamical model calculations predict that V 2 as 
a function of pt follows mass ordering, where the V 2 of 
heavier hadrons is lower than that of lighter hadrons and 
vice-versa [3, [H, [13 . Mass ordering is indeed observed 
in the identified hadron measured in the low px 
region {px <1.5 GeV/c) Recent phenomenological 
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FIG. 4. (Color online) V 2 {(j))/v 2 {p) ratio as function of pt for 
0-30%, 30-80% and 0-80% centrality in Au + Au collisions 
at -^snn = 200 GeV. Shaded boxes are the systematic un¬ 
certainties and vertical lines are the statistical uncertainties. 
The first data point of 0-80% centrality is shifted towards 
right by 400 MeV/c. The bands in panel (a) and (b) repre¬ 
sent the hydro model |3 and transport model calculations for 
V 2 { 4 >)/v 2 {p), respectively. 


calculations, based on ideal hydrodynamics together 
with a hadron cascade (JAM), show that the mass 
ordering of V 2 could be broken between cj) mesons and 
protons at low px {px <1.5 GeV/c) [3,13 ■ The broken 
mass ordering is thought to be due to late stage hadronic 
re-scattering effects on the proton V 2 , since the model 
calculations assume a low hadronic cross-section for the 
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(j) but a large hadronic cross-section for the proton. 

The ratios of (j) V 2 and proton V 2 are shown in Figured 
The ratios are larger than unity at pr ^ 0.5 GeV/c for 
0-30% centrality showing an indication of breakdown of 
the expected mass ordering in that momentum range. 
This could be due to a large effect of hadronic re¬ 
scattering on the proton V 2 , indicated by the shaded red 
band in panel (a) of Figure 01 We have also considered 
the effect of the momentum resolution and energy loss 
of the TPC as well as decay (feed-down) effects on the 
proton V 2 - Our study, based on the UrQMD framework, 
indicates that the momentum resolution and decay 
effects on the ratio of V 2 { 4 >) to V 2 {p) in the measured 
momentum region are negligible. The break down of 
mass ordering of V 2 is more pronounced in 0-30% than 
in 30-80% centrality. For example, the ratio U 2 (((')/u 2 (p) 
is 4.35 ± 0.98±!];t| at pr = 0.52 GeV/c in 0-30%, while 
it is 1.12 ± 0.10 ±q;q|3 in 30-80%. In the central events, 
both hadronic and partonic interactions are larger than 
in peripheral events. Therefore, the combined effects of 
large partonic collectivity on the (j) V 2 and significant late 
stage hadronic interactions on the proton V 2 produce 
a greater breakdown of mass ordering in the 0-30% 
centrality data than in the 30-80% [l^. This observa¬ 
tion indirectly supports the idea of a small hadronic 
interaction cross-section for the 4> meson. We have also 
studied the ratio of 4> V 2 to proton V 2 using the transport 
models AMPT and UrQMD [^. The U 2 ((/))/u 2 (p) 
ratio for 0-30% centrality from AMPT and UrQMD 
model are shown in Figure 0] (panel b). The black 
shaded band is from AMPT with a hadronic cascade 
time of 0.6 fm/c while the yellow band is for a hadronic 
cascade time of 30 fm/c. It is clear from Figure 0] (panel 
b) that with increasing hadronic cascade time (and 
therefore more hadronic re-scattering), the U 2 ((/))/u 2 (p) 
ratio increases. This is attributed to a decrease in the 
proton V 2 due to an increase in hadronic re-scattering 


while the meson V 2 remains unaffected [l5j| . The ratios 
from UrQMD are shown as a red shaded band which 
is much smaller than unity. The UrQMD model lacks 
partonic collectivity and therefore does not fully develop 
the (/-meson V 2 - 

In summary, we have reported high-statistics elliptic 
flow measurements for multi-strange hadrons (S and VL) 
and (j) meson with other light and strange hadrons (tt, AT, 
, p and A) in Au -|- Au collisions at ^snn = 200 GeV 
for different centralities. The pt dependence of (j) and 
V 2 is observed to be similar to that of tt and p, indicating 
that a large amount of collectivity is developed in the 
initial partonic phase for light and strange hadrons. 
NGQ scaling holds within the statistical uncertainty 
for both 0-30% and 30-80% centralities, suggesting 
collective motion of quarks prior to hadronization. The 
comparison between the (j) and p V 2 shows that at low pr , 
there is a possible violation of hydrodynamics-inspired 
mass ordering between (j) and p. Model calculations 
suggest that the pr dependence of V 2 { 4 >)/V 2 (jp) can 
be qualitatively explained by the effect of late-stage 
hadronic re-scattering on the proton V 2 0:0 ■ 
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